Highly polarized Fermi gases in two dimensions 



Meera M. Parish^' ^ and Jesper Levinsen^ 

^London Centre for Nanotechnology, Gordon Street, London, WCIH OAH, United Kingdom 
^Cavendish Laboratory, J J Thomson Avenue, Cambridge, CBS OHE, United Kingdom 

(Dated: January 17, 2013) 

We investigate the highly polarized limit of a two-dimensional (2D) Fermi gas, where we effectively 
have a single spin-down impurity atom immersed in a spin-up Fermi sea. By constructing variational 
wave functions for the impurity, we map out the ground state phase diagram as a function of mass 
ratio m^/mi and interaction strength. In particular, we determine when it is favorable for the 
dressed impurity (polaron) to bind particles from the Fermi sea to form a dimer, trimer or even 
larger clusters. Similarly to 3D, we find that the Fermi sea favors the trimer state so that it exists 
for m^/mi less than the critical mass ratio for trimer formation in the vacuum. We also find a 
region where dimers have finite momentum in the ground state, a scenario which corresponds to 
the Fulde-Ferrell-Larkin-Ovchinnikov superfiuid state in the limit of large spin imbalance. For equal 
masses (m^ = ^i)^ we compute rigorous bounds on the polaron-dimer transition, and we show 
that the polaron energy and residue is well captured by the variational approach, with the former 
quantity being in good agreement with experiment. When there is a finite density of impurities, we 
find that this polaron-dimer transition is preempted by a first-order superfluid-normal transition at 
zero temperature, but it remains an open question what happens at finite temperature. 



I. INTRODUCTION 

The creation of tunable Fermi systems with ultracold 
atomic gases has greatly renewed interest in fermionic 
pairing phenomena and superfluidity. In particular, re- 
cent experiments have successfully confined fermionic 
atoms to a quasi- two- dimensional (quasi-2D) geome- 
try [1-7], thus enabling the investigation of 2D Fermi 
gases. Such model 2D systems can potentially provide 
insight into more complicated solid state systems such as 
electron-hole bilayers [8, 9] and high temperature super- 
conductors [10]. Moreover, 2D is the "marginal" dimen- 
sion where quantum fluctuations are enhanced compared 
to 3D, leading to the destruction of Bose Einstein conden- 
sation (BEC) at finite temperature. At the same time, 
there is a dearth of exactly solvable models, in contrast 
to the case in ID. This makes the cold-atom system even 
more important for testing theoretical approaches in 2D. 

A canonical many-body problem currently under inves- 
tigation is that of an impurity interacting with a medium. 
This so-called "polaron" problem traditionally involves a 
bosonic medium such as a bath of phonon excitations in 
a crystal [11]. However, the advent of cold atoms has 
extended this problem to include a fermionic medium, 
where the impurity can now undergo a sharp transition 
and effectively change its statistics by binding fermions 
from the background Fermi gas [12, 13]. For a fermionic 
impurity, this corresponds to the highly polarized limit 
of a two-component (t, i) Fermi gas and, thus, the exis- 
tence of such binding transitions will determine the topol- 
ogy of the phase diagram for the spin-imbalanced Fermi 
gas [14, 15]. In 3D, the | impurity is initially dressed 
by density fluctuations of the t Fermi gas, forming a 
polaron state [16], but with increasing attractive inter- 
actions it can bind an extra t particle to form a dimer 
(molecule) [12, 13, 17-20] or it can even bind two t par- 
ticles to form a trimer [20], depending on the mass ratio 



m^/m^. For the 2D case, there has been debate about 
whether strong quantum fiuctuctions preclude the exis- 
tence of such binding transitions [21, 22], but one of us 
has recently argued that sharp polaron-molecule tran- 
sitions can exist [23] and this appears consistent with 
experiment [6, 24]. 

In this work, we further investigate this issue by map- 
ping out the ground state phase diagram for the 2D im- 
purity problem as a function of mass ratio m^/m^ and 
interaction strength. Similarly to 3D, we show that the 
impurity can bind one, two or more t particles from 
the Fermi sea to form bound clusters dressed by den- 
sity fluctuations. We approximate the different impu- 
rity states using variational wave functions that include 
a finite number of particle-hole excitations of the Fermi 
sea. We find that the presence of a Fermi sea appears 
conducive to the formation of trimers, which in vac- 
uum exist for m^/m^ > 3.33 [25]. On the other hand, 
tetramers containing three t particles and existing for 
m^/mi > 5.0 in vacuum [26], appear to be disfavored 
in the many-body system. We further find a region of 
parameter space in which the impurity binds a single 
particle from the Fermi sea to form a dimer with finite 
momentum. A finite density of such molecules will form a 
spatially modulated condensate and thus this state is the 
single-particle analog of the elusive Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) phase [27, 28]. 

The presence of bound clusters of + 1 particles such 
as trimers and tetramers (A^ = 2 and 3) introduce impor- 
tant (A^ + l)-body correlations. Even if these states are 
metastable, their presence can still affect the system: for 
instance, we expect three-body correlations to be more 
important in 2D than in 3D for the mass-balanced system 
due to the smaller repulsive barrier between the iden- 
tical fermions, and the resulting enhanced three-body 
interactions [29, 30]. In order to quantify the effect of 
these higher order correlations, we consider the polaron 
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state dressed by two particle-hole pairs for equal masses 
(m^ = m^). This wavefunction includes three-body cor- 
relations, whereas the usual Chevy wavefunction [16] con- 
tains only two-body correlations. The results of these two 
variational wavefunctions match in the limit of weak in- 
teractions and we also find a reasonable agreement across 
the regime of strong interactions. This approach further 
allows us to arrive at rigorous upper and lower bounds 
for the position of the polaron-molecule transition, yield- 
ing -0.97 < ln(A:Ft^i2D) < -0.802, with kp^ the Fermi 
momentum of majority t particles and a2D the 2D scat- 
tering length. 

At any given interaction strength, only one state may 
be the ground state and all other states are at best 
metastable. The polaron wavefunction of Ref. [16] leads 
to a completely real equation for the polaron's energy, 
and thus one might think that the variational approach 
may only be used to calculate ground state properties 
of the polarized Fermi gas. On the other hand, using 
a diagrammatic approach as in Ref. [31] to describe the 
polaron, naturally incorporates an imaginary part in the 
energy, and thus the metastability of excited states is 
automatically present. We resolve this seeming discrep- 
ancy by showing here that a straightforward change of 
the variational wavefunction and the minimization pro- 
cedure allows the variational approach to be extended 
to describe not only the ground state but also excited 
branches. Thus the two approaches are, in fact, fully 
equivalent. 

The experimental observation of the transition from a 
polaron to a molecule may be precluded by phase sepa- 
ration into an unpolarized superfluid phase and a fully 
polarized normal phase. Using results from a quantum 
Monte Carlo (QMC) simulation of the superfiuid phase in 
an equal-mass unpolarized Fermi gas [32] , we indeed find 
this to be the case in the highly-polarized limit. How- 
ever, this result only strictly applies to zero temperature 
and it remains an open question whether this occurs at 
finite temperature. 

The paper is organized as follows: The contact inter- 
action is introduced for the 2D system in Sec. II. Sec- 
tions III and IV introduce the variational wave functions 
used to describe the possible ground states of the impu- 
rity atom, while Sec. V describes the resulting single- 
impurity phase diagrams. In Sec. VII we derive the 
condition for phase separation precluding the polaron- 
molecule transition, while in Sec. VI we introduce the 
polaron dressed by two particle-hole pairs. In Sec. VIII 
we show the manner in which the variational method may 
be modified to describe metastable states, and finally in 
Sec. IX we conclude. 



II. PRELIMINARIES 

We consider the Hamiltonian for a two-component 
atomic Fermi gas interacting via a short-range interac- 



tion in two dimensions (2D): 



k,k',q 



where the spin a =t, |, eka 



2mcT 



, and g is the strength 



of an attractive contact interaction. We work in units 
where h and the system area are both 1. Note that since 
we are considering low-energy, 5- wave interactions, the 
Pauli exclusion suppresses interactions between the same 
species of fermion. For the two-body problem (t and 1) , 
we simply have 
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g £b + Ckf + £14 



(2) 



where A is the UV cut-off and 5^ is the binding energy 
of the weakly bound diatomic molecule which always ex- 
ists for an attractive interaction in 2D. We see that the 
integral logarithmically diverges if we fix sb and take 
A ^ 00. Thus A cannot be removed from the problem 
and the binding energy sb depends on both A and in 
contrast to ID. For the many-body system, our results 
become independent of the cut-off A once Eq. (2) is used 
to replace g with sb- 

At low energies, the elastic scattering of a spin-t and a 
spin- 4. atom at relative momentum q is described by the 
s-wave scattering amplitude [33] 



27r 



ln[l/(ga2D)] +i7r/2' 



(3) 



Here, the momentum dependence of the scattering ampli- 
tude is characterized by the 2D scattering length a2D > 0. 
Analytic continuation of the scattering amplitude to mo- 
menta on the positive imaginary axis yields a pole in the 
scattering amplitude at g = ia2^. This corresponds to 
the binding energy Sb = l/2mr-a2i3, with the reduced 
mass defined as rrir = m^mi/{m^ -h Ta^). Thus, the 2D 
scattering length essentially corresponds to the size of 
the two-body bound state. Unlike in 3D, the scattering 
amplitude in 2D does not reduce to a constant at low 
scattering momenta; instead the scattering is strong in 
the regime ln(g'a2D) ~ and weak when | ln(ga2D)| ^ 1- 
In the presence of a Fermi sea, the characteristic momen- 
tum scale is set by the Fermi momentum kp and we thus 
expect strong many-body effects when ln(/ci?<^2D) ^ 0. 
Indeed, Bloom [34] has demonstrated how the 2D Fermi 
gas in the regime | ln(A:i?a2D)| ^ 1 is perturbative in 
l/ln(A;Fa2D). 

In the following, we focus on the case where we have 
a single spin-down minority atom immersed in a non- 
interacting Fermi sea of spin-up atoms — the extreme 
limit of population imbalance. Note that the phase di- 
agram we obtain for the single impurity atom is inde- 
pendent of impurity statistics and is thus also relevant 
to 2D Bose- Fermi mixtures. However, the focus of this 
paper will be on fermionic impurities. Defining the inter- 
action parameter in the imbalanced gas, 77 = \n{kF^a2-D), 
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in the limit of weak attractive interactions (or, equiva- 
lently, the large- density limit) where ^ 1, the ground 
state is expected to approach that of the non-interacting 
system: |Pi(0)) = Cq^ I^^*), where \FS) represents the 
Fermi sea of t-particles. Note that the subscript on the 
"polaron" wave function \P) denotes the number of op- 
erators acting on the Fermi sea; this is the nomenclature 
we will use throughout the paper. Decreasing ry will even- 
tually give rise to one or more binding transitions, where 
the spin-down minority particle binds one or more spin- 
up particles. The behaviour of the system in this regime 
can be analysed with the use of variational wave func- 
tions for the different bound states. 



III. "UNDRESSED" WAVE FUNCTIONS 

In this section, we neglect the particle-hole excitations 
generated by the impurity ^.-particle interacting with the 
Fermi sea and focus on the "bare" part of the exact wave 
function. For the molecule and polaron, this is equivalent 
to the mean field approach for the spin-imbalanced Fermi 
gas (see, e.g., Ref. [35]). In the present problem of a sin- 
gle impurity interacting with a Fermi gas, the approach 
serves as a useful introduction to the more complicated 
wavefunctions dressed by particle-hole fluctuations. It 
should also provide insight into the relevant few-body 
correlations present in the system. Below, in Sec. IV, 
we consider wave functions dressed by one particle-hole 
excitation, an approach which yields quantitatively ac- 
curate results in the perturbative regimes of weak and 
strong attractive interactions and which is also expected 
to give a good description across the regime of strong 2D 
interactions 77 ~ 0, as we discuss below. 

The results of this section become exact in the vac- 
uum limit, kp'i- ^ 0. In this case, we can study tran- 
sitions between different few-body states containing a 
single spin-^ atom and spin-^ particles. In vacuum, 
a bound diatomic molecule always exists for two atoms 
interacting via contact interactions under a strong two- 
dimensional conflnement. Additionally, a trimer state 
consisting of two heavy fermions and one light par- 
ticle becomes energetically favorable for a mass ratio 



r = m^/m^ above 3.33 [25], while a tetramer containing 
three heavy fermions and one light particle is favorable 
for r > 5.0 [26]. 

A. Molecules (M2) 

The simplest, lowest-order variational wave function 
for a bound pair or "molecule" is: 

|M2(p)) = ^V'L^)c;_k,4,|F5), (4) 

k 

where p corresponds to the center-of mass momentum of 
the pair, while the spin- up particle momentum satisfies 
k = |k| > kp^. This wave function gives the exact two- 
body state in the limit /ci?^ ^ and is in fact identical to 
the BCS mean field wave function for extreme imbalance 
(after the state with one | particle has been projected 
out). Note that we assume the wave function (4) has one 
less t particle in the Fermi sea compared with the non- 
interacting state I Pi) in order to preserve particle num- 
ber. Indeed, we will assume throughout this paper that 
all the impurity wave functions we introduce have the 
same number of t particles as |Pi ). This is equivalent to 
measuring the energy of t particles with respect to their 
chemical potential, the Fermi energy Ep^ = k'^^/2m^^ 
and thus we define ^kt = ^kt — Ep^. 

Minimizing the expectation value (M2(p)| {H — 
E) |M2(p)) with respect to ip^^ yields an implicit equa- 
tion for the molecule energy E: 



9 ^-^ + ekt + ep_k;' 

Here and in what follows we use the convention that the 
momentum k corresponds to a particle excited out of the 
Fermi sea, i.e. \k\ > kF'[- Additionally, the energy E is 
assumed to be with respect to the (macroscopic) energy 
of the non-interacting t Fermi gas. We neglect Hartree 
terms involving gkp^/4:7r since these vanish when we take 
the limit A — > 00, ^ ^ 0. Converting sums into integrals 
and sending A ^ 00 then gives: 



2sB = -E' 



p^{r — l)r 



2m^(l-hr)2 l-hr\ 



kl 



Ft (p -kp^y 



2mr 



2mr 




E'^lE'^ 



2mr 



2mr 



where we defined the energy E' = E — Clearly, 

when we approach the two-body limit (/ci?^ ^ 0)? the 
molecule has its lowest energy at zero momentum. For 

this case, we simply have E = 2^ — Sb- However, once 



^Ft^2D > l/v^5 we find that the molecule acquires a 
finite momentum: 

1 -h r I 

P = V/^Ft«2D\/r - 1- (6) 

Increasing /ci?t^^2D further eventually causes the molecule 
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to unbind into the state |Pi(0) ). At this unbinding tran- 
sition (assuming there is a direct transition), we always 
find that p = kp^, in contrast to the 3D case where p = 
for mass ratios sufficiently close to one [20]. Referring 
to Eq. (6), this means that the molecule unbinds when 
^Ft^2D = (1 + i-G- when sb equals the center-of- 

mass kinetic energy of the molecule ai p = kp^. 



B. Trimers (Ts) 

Another possible bound state is the trimer consisting of 
two spin-up fermions and one spin-down particle. Since 
this trimer involves identical fermions, its angular mo- 
mentum L must necessarily be odd and so the lowest- 
energy trimer is ap-wave (L = 1) bound state. Indeed, it 
may be regarded as a p-wave pairing of spin- up fermions 
mediated by their 5- wave interactions with the spin-down 
particle. The lowest-order variational wave function for 
the trimer is 



1^3(0))= ^7k,k.cl 

kik2 



ki 



(7) 



Since angular and linear momentum do not commute, 
we restrict ourselves to trimers with zero center-of-mass 
momentum (p = 0) in order to have a well-defined L. 
Such a restriction is unlikely to be drastic since there is 
no physical reason to believe that the trimer will have 
its lowest energy at finite momentum. In fact, we would 
generally expect the energy to be higher at non-zero p 
since the wave function would then contain an admixture 
of higher-energy angular momentum states L > 1. 

Minimizing {T^l {H — E) \Ts) and defining the function 
/k2 = Ski ^kik2 then results in the equation 



yk2 



E 



^kik2 



E 



/ki 
^kik2 



(8) 



where Ek^ks = + ^ki+ksi + ^kit + ^kst- For L = 1, 
we have /k = fk^^^i where (j) is the angle with respect 
to the X-axis and / is an arbitrary function. Performing 
the angular integration, leaves a one-dimensional integral 
equation for fk which is subsequently solved by discretiz- 
ing /c-space and converting the integral equation into a 
matrix eigenvalue equation [36]. We emphasize that at 
this order of approximation, the effect of the Fermi sea is 
only taken into account in Eq. (8) through the restriction 
on the momenta in the sums (|ki| > kp^)- The energy 
of the trimer in the absence of a Fermi sea is recovered 
upon taking the limit kp'i ^ 0. 



spin-t particles: 

l^^+l)= Xki...k^cLk,..._k^^4,f ••^Ltl^^)• 

kl...k^^ 

(9) 

Again, minimizing (Xtv+iI {H — E) \Xn+i) and defining 
/ks-.k^r = Eki Xki...kAr yields an equation for the energy 
of the -h 1 bound state: 



/k2...i< 



^ E^ 



ki...kA 



/kik3...kAr H- /k2kik4...kjv + • • • + /k2...kiv_iki 



ki 



£^1 



ki...kjv 



where £;ki...kiv = -£ + eki+...+kAr; + ^k^f This equa- 
tion is explained in detail in Appendix A where an alter- 
native derivation in terms of diagrams is presented. 

Again the vacuum limit is recovered by letting kp^ ^ ^ 
and we stress that Eq. (10) is quite general in this limit: 
It is equally valid for ID, 2D and 3D systems; it may be 
extended to narrow Feshbach resonances by letting the 
coupling constant be energy dependent; and it may be 
used to treat the quasi-2D problem as in Ref. [26], if one 
includes a summation over harmonic oscillator modes. 
The equation satisfied by the tetramer energy {N = 3) 
in this limit was obtained for the 3D problem in Ref. [37] , 
the equation for the A/" + 1 problem in a quasi-2D geom- 
etry was derived in Ref. [26]. Finally, Ref. [38] derived 
an expression similar to our Eq. (10) for the 3D N -\- 1 
vacuum problem [39] and Ref. [40] generalized this result 
to include also ID and 2D. 



IV. WAVE FUNCTIONS WITH ONE 
PARTICLE-HOLE EXCITATION 



We can improve on the wave functions in Sec. Ill by 
adding a single particle-hole excitation on top of the 
Fermi sea. We expect these improved wave functions 
(which are perturbative in the number of particle-hole 
excitations) to provide a reasonably accurate estimate of 
the single- impurity energy even in the regime of strong 
interactions, kp^a^D ^ 1, since it has been argued that 
contributions from two or more particle-hole excitations 
nearly cancel out via destructive interference [41]. We 
further investigate the validity of this approach in Sec. VI 
where we present our results for the impurity dressed by 
two particle hole pairs. 



C. The N+1 problem 



A. Polarons (P3) 



The above approach may obviously be extended to the 
bound state consisting of a single spin-| impurity and N 



Adding one particle-hole excitation to the non- 
interacting state I Pi ) gives the improved polaron wave 
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function: 

mp)) = \FS) + E"S4+q-k,4tCqt i^^^), 

kq 

(11) 

where we have now included a center-of-mass momentum 
p. Minimizing {Ps{p)\{H — E)\Ps{p)) then gives us 
equations: 



kq 

,.(P) _ 



(12) 



{E - ep+q_k; - ekf + Cqf ) a^q - 

5E"S-5E<'+^^"o''^ (13) 



We remind the reader that particle momenta k satisfy 
|k| > kF^- The hole momenta will be denoted by q with 
|q| < kp^. Combining Eqs. (12) and (13) gives an im- 
plicit equation for the energy, 



1 V- 1 



(14) 



with £;pqk = -£^ + ^kt-^qt + ^P+q-H- Eq. (14) has two 
solutions: The attractive and repulsive polaron which 
have energies above and below the energy of the free im- 
purity, respectively. Note that we must include here by 
hand a small imaginary part on the right hand side in 
order to describe the metastable repulsive polaron. We 
show in Sec. VIII how to the extend the variational ap- 
proach so that such imaginary terms arise naturally for 
metastable states. 

The solution of Eq. (14) at small but finite momentum 
|p| <C kp^ gives the dispersion: 



^(P) = ^(0) 



2m* ' 



(15) 



with m* the effective mass. In Fig. 1 we display the 
effective mass of the attractive polaron for three different 
mass ratios. Here, r = 6.64 (r = 1/6.64) corresponds to 
the experimentally relevant lithium (potassium) impurity 
in a potassium (lithium) Fermi sea. As opposed to 3D, 
where the effective mass was found to diverge once the 
polaron was a metastable excitation [17], we find that the 
effective mass is always positive in 2D, i.e. the polaron 
always has its minimum energy at p = 0. 

In addition to the effective mass, the polarons are also 
described by the wave function overlap with the free im- 
purity — the residue: 



Z = a, 



(16) 



The residue of the attractive polaron goes to 1 in the limit 
^Ft<^2D ^ 1 where the polaron is a well defined quasipar- 
ticle, while it vanishes in the opposite limit. The repul- 
sive polaron displays the opposite behavior. We discuss 





1 1 1 1 1 1 1 1 1 1 1 


1 1 1 














- 




- 










\ 






\ 




■ 


\ 






\ 






\ 




_ 


\ 


_ 




\ 






\ 




- 


\ 


- 



















































FIG. 1: (Color online) Effective mass of the attractive po- 
laron calculated within the variational ansatz P3 as a func- 
tion of interaction 77 = ln(A:FtQ^2D) for different mass ratios: 
m^/mi=6M (dashed), 1 (solid), 1/6.64 (dotted). For the 
equal mass case, our results match those obtained in Ref. [42]. 
The ground state is a polaron for weak attractive interactions 
but becomes a molecule (trimer) to the left of the diamond 
(triangle). Here we compare the wavefunction P3 with M4 
and T5, Eqs. (17) and (21). 



the residue of the attractive polaron further in Sec. VI 
where we compare with the result of dressing the impu- 
rity by two particle-hole pair excitations. 

The polaron state has been thoroughly studied in 
the 2D geometry. In particular, theoretical studies of 
IP3) have obtained the energy of the attractive polaron 
[21, 23]; the effective mass and residue of the attractive 
and repulsive polarons [42, 43]; additionally, Ref. [43] 
found the lowest lying repulsive polaron state in the limit 
kFt^2D ^ 1 to have finite momentum. Experimental ev- 
idence of both the repulsive and the attractive polaron 
was found using radiofrequency spectroscopy in Ref. [6]. 



B. Molecules (M4) 



We can likewise improve on the bare molecule state 
IM2) by adding one particle- hole pair as follows: 



iM4(p)) = E^lrH-ki4ti^5) 

k 

+ E ^uu'A+^-^-^'AA't^^t \FS). (17) 

kk'q 
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FIG. 2: (Color online) Inverse effective mass of the molecule 
as a function of interaction parameter for different mass ratios: 
mt/m^= 1 (solid), 1/6.64 (dotted). 




FIG. 3: Illustration of the dispersion of M^ip) for = Sm^.. 
(a) Momentum at the minimum of the dispersion as a function 
of interaction parameter, (b) Dispersion at 77 = —1.25 and (c) 
dispersion at 77 = —1.1. 



The minimization procedure now leads to the following 
equations for the energy of the molecule: 



^pqkik 



E 



1 

^qk 



E 

kq 



^qk 



(18) 



pqkik 



pki 



E 



Gaf 



q'ki 



pki 



(19) 



where E^i^ = -E ^ + ^p+k; and £^pqkik2 = -E ^ 
^kit + ^^k2t - W + ep+q_ki-k2i- The function G is de- 
fined as Gqk = Ek' ^kk^'q/Ek' ^k^^- The energy of 
the molecule in the 2D geometry was first obtained in 
Ref. [23]. The coupled integral equations (18) and (19) 
are equivalent to the equation for the molecule energy 
obtained in Refs. [17-19] for an impurity in a 3D Fermi 
gas. 

Similarly to the polaron above, the energy of the 
molecule as a function of momentum yields the dispersion 



Eip) = EiO) 



P 



2M* 



(20) 



with M* the effective mass. In the absence of inter- 
actions, the effective mass of the molecule is simply 
M* = + m^. However, we observe a strong depen- 
dence of the effective mass on the interaction parameter, 
as shown in Fig. 2. As further illustrated in Fig. 3, once 
M* < the molecule has its minimum energy at finite 
momentum, i.e. the pairing occurs at a finite momen- 
tum. As mentioned previously, this corresponds to the 
FFLO phase in the limit of large polarization, and we 
find that this phase occupies regions of the phase dia- 
grams discussed in Sec. V. 



C. Trimers (T5) 

As discussed above, in the absence of a Fermi sea and 
for r > 3.33, the impurity binds two particles to form a 
trimer. Since we wish to investigate the possibility of the 
trimer being the ground state across the regime of strong 
interactions, we use the dressed wave function 



kik2 



E 

kik2kq 



7kik2kqC^_k,-k2-ki4it^k2t^kt^qt 1^^)- (21) 



Following the minimization procedure, we find two cou- 
pled integral equations: 



ki 



1 



^kik2 



E 

k2 



^k2 
^kik2 



qkik2 
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■E 



E' 



1 



^ £^qkik2k 
Gqkik + Gqkk 



E 

qk2 



qkik2 
^kik2 

-^k2 



(22) 



E. 



qkik2k 



E 

q' 



^kik2 
Gq'kik2 



^kik2 



(23) 



We have defined Jk = ^Ek'7kk' and Gqkik2 = 
3^Ek' 7kik2k'q. The energies are ^kik2 = -E ^ ^kit + 
^^k2t + ^ki+k2i and ^qkik2k = -E + ^kit + ^k2t + ^kt - 
^q^ + eq_k-ki-k2i- Eqs. (22) and (23) were first derived 
and solved for the trimer energy in Ref. [20] for the im- 
purity problem in a 3D Fermi gas. The projection onto 
the p-wave trimer state is performed by taking 

/ki = 4e*^S Gqk.k. = G{q,h,k2,A4>i„A4>2g)e'^\ 

(24) 

with ^1, 02, and (j)q the angles which ki, k2, and q make 
with the axis of reference, while A^ig = — (j)q and 

A02q = (/>2 - (/>g. 
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FIG. 4: (Color online) Ground state phase diagram for an 
impurity atom attractively interacting with a 2D Fermi gas. 
The FFLO phase corresponds to the molecule M2(p) with 
non-zero momentum p in the ground state. The behavior of 
the momentum at which the energy is at its minimum is given 
by Eq. (6) and it goes smoothly to zero at the dashed (red) 
line given by m^/mi = l/(A:Fa2D)^. We find that small slivers 
of FFLO and trimer phases remain as 77 ^ oc. Note, also, 
that the trimer exists above the critical mass ratio r ^ 3.33 
in the limit 77 ^ —00, which agrees with the result for the 
3-body bound state in a vacuum [25]. 



D. The AT + 1 problem 

Like the bare wavefunctions of Sec. Ill, the above ap- 
proach may be extended to the study of the bound states 
of the impurity and N spin-^ fermions. The variational 
wave function is dressed by one particle-hole pair exci- 
tation of the Fermi sea and the minimization procedure 
carried out as above. This leads to two coupled integral 
equations similar to Eqs. (22) and (23) for the trimer 
above. The equations are derived in Appendix B using 
the diagrammatic technique. We shall not attempt here 
to solve for the energy of the tetramer or bound states 
containing even more particles. 



V. PHASE DIAGRAMS 

We now determine the ground state for the single 
impurity and the correponding binding transitions. In 
Fig. 4 we show the phase diagram for the "undressed" 
wave functions of Sec. III. Surprisingly, we find that a 
molecule existing at a given mass ratio r < 3.33 must al- 
ways first bind an extra spin-up fermion to form a trimer 
before it can unbind into a polaron. This appears to be 
an artifact of the approximation. However, it does signify 
the importance of three-body correlations for all mass ra- 
tios in 2D. Additionally, we find a sliver of FFLO phase, 
corresponding to a finite momentum molecule, on the 
border of the zero-momentum molecule and the trimer 



FIG. 5: (Color online) Ground state phase diagram for a 
spin-down impurity atom immersed in a 2D Fermi gas, with 
phase boundaries calculated using the dressed wave functions 
of Sec. IV. The FFLO phase corresponds to the molecule 
M4(p) with non-zero momentum p in the ground state. See 
Fig. 3 for the behavior of p as a function of 77 across the FFLO 
region. 

phases. The large region of trimer phase below r 3.33 
appears to result from the fact that the FFLO molecule 
is unstable towards binding an extra t particle, like in 
3D [20]. Whereas trimers are favored by the medium, 
we find that tetramers consisting of three spin-^ parti- 
cles and the impurity appear to be disfavored, i.e. the 
phase transition is found to occur at larger mass ratios 
in the medium than the critical mass ratio of r = 5.0 in 
vacuum [24]. This suggests that four-body correlations 
are not as important in the many-body system at low 
mass ratios as might be initially expected. 

Next, Fig. 5 shows our phase diagram obtained using 
wave functions dressed by one particle-hole pair - see 
Sec. IV. As above, we find that the trimer is favored 
by the Fermi sea, but it now does not appear below a 
mass ratio of r ~ 2.1. Additionally, we find that the 
FFLO region, where the ground-state molecule has finite 
momentum, is enlarged at this level of approximation, 
making it possible that the FFLO phase may be observed 
in this system. 

We expect this phase diagram to be qualitatively cor- 
rect also across the regime of strong many-body correc- 
tions, |1/?7| <C 1, as contributions from two or more 
particle-hole pairs cancel approximately [41] and parts 
of the phase diagram are fixed by perturbative and exact 
calculations. In the limit of large negative 77, the dress- 
ing of the molecule and trimer by one particle-hole pair 
yields the correct form of the first-order correction to 
their energy due to the interaction with the Fermi sea in 
a perturbative expansion in 1/|77|. Thus the phase tran- 
sition from M4(0) to T^{0) and the fact that the trimer 
is favored by the Fermi sea is quantitatively robust. Fur- 
thermore, for an impurity with infinite mass (r = 0), 
the approach correctly identifies the dimer as the ground 
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FIG. 6: (Color online) Energy measured from the two-body 
binding energy. The dashed line corresponds to the results for 
the ansatz IPs), Eq. (11), the solid hne to IP5 ), Eq. (25), and 
the dot-dashed line to the molecule ansatz IM4), Eq. (17). 
The polaron-molecule transitions in the two approximations 
are illustrated by vertical dotted lines, while the experimental 
data is taken from Ref. [6], rescaled to reflect the quasi-2D 
nature of the experiment [44]. 



state [23]. 



VI. ACCURACY OF THE VARIATIONAL 
APPROACH 

The variational approach employed above has been val- 
idated by several methods in various settings. In 3D, the 
energy and effective mass of the attractive polaron quasi- 
particle have been measured in Refs. [45] and [46], respec- 
tively, and good agreement with the variational calcula- 
tion [16, 31] was obtained, even close to the unitary limit. 
It has also been shown that the polaron and molecule 
energies [17-19, 41] calculated using this method for a 
3D Fermi gas with = are in good agreement with 
those from quantum Monte Carlo [12, 13]. More recently, 
an experimental and theoretical investigation of the pop- 
ulation imbalanced ^^K-^Li mixture showed impressive 
agreement between the measured energies, residues, and 
lifetimes of the attractive and repulsive polarons when 
compared with the theoretical predictions from the dia- 
grammatic technique [47] (which is equivalent to the vari- 
ational one presented here). Finally, it has been argued 
that contributions from two or more particle-hole excita- 
tions nearly cancel out via destructive interference [41]. 
This argument does not depend on the dimension and in- 
deed in the one-dimensional geometry it was shown that 
the variational ansatz agrees well with the exact Bethe 
Ansatz solution [48]. 

We now wish to demonstrate explicitly the accuracy of 
the dressed wave functions of the previous section. We 
emphasize the perturbative nature of the many-body sys- 
tem as long as 1/| ln(/ci?^a2D)| ^ 1, but we wish here 



to quantify the effect of quantum fluctuations in the 
strongly interacting region, kF^CL2B ~ 1. To this end we 
write down a variational wave function for the impurity 
dressed by two particle-hole pairs: 

|P5(0) ) = ao4 \FS) + Y, akq<_k^4^cqt \FS) 

kq 



E 

kik2 



<^kik2qiq2<^qi+q2-ki-k2i^kit^k2t^^i^^^2t 1^^) 



(25) 

This wave function was first studied for the 3D polaron 
problem in Ref. [41]. The minimization of (P5I {H — 
E) IP5) then results in two coupled integral equations, 



1 

qk 



q' 



kq' 



qk 



(26) 



qiq2k 



^qiq2kk' 



E 



qiq2k' 

^qiq2kk' 



/qi + Xlq'2 Gqiqsk /q2 Zlq^ Gq/^q2k 



E, 



qik 



E. 



q2k 



(27) 



The energies are £^qk = —E^ ~ W + ^q-k^ and 

^qiq2kik2 = " ^ + ^kit^'k2t ~ ^'qit ~ ^q2t + ^qi+q2-ki -k2i 7 

and we define Gq,q2k = 4^ Z)k' <^k'kqiq2 and /q = 
g Qfkq- The restriction to s-wave scattering of the im- 
purity off a majority atom means that /q depends only 
on the magnitude of q. Likewise, the vertex G depends 
on the magnitudes of qi, q2, and k and the two angles 
^qi,k and Zq2,k. 

Our results for the polaron energy for equal masses 
of the two species are shown in Fig. 6. We see that 
the two ansatze agree very well in the weakly interact- 
ing regime, 77 ^ 1. In the regime of strong many-body 
effects, 1/1 7^ I <C 1, the polaron energies resulting from 
the two ansatze show the same non-monotonic behavior 
and are never further separated than by {).1Ef^. The 
polaron in a quasi-2D geometry was observed in a recent 
experiment [6] . The data from the experiment is also dis- 
played in Fig. 6 and is seen to match the energy of the 
IP5 ) ansatz quite well for 7^ > 0.25, leading us to believe 
that the ) ansatz is likely to closely reproduce the true 
energy of the polaron. The reason for the discrepancy be- 
low this value may be due to the finite interaction range 
of the interatomic potential, a finite density of impurity 
atoms, finite temperature effects, and trap-averaging ef- 
fects. 

From the polaron energies in the two ansatze we obtain 
the following values of the interaction parameter at the 
polaron-molecule transition: 



P3-M4 transition: rj 
P^-M4 transition: 77 



-0.80, es/EFt 
-0.97, es/EFt 



9.9 
14.0 



(28) 



This compares well with the experimental result of 
Koschorreck et al [49], r] = -0.88(0.20) [44]. 
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If the polaron and molecule variational wavefunctions 
converge with increasing numbers of particle-hole pair ex- 
citations, then it is possible to show that the above results 
for the polaron-molecule transition, Eq. (28), provide up- 
per and lower bounds for the actual polaron-molecule 
transition. First, we define ^2n-i — and — 
to be the energies for the polaron \P2n-1) and molecule 
|M2n ), respectively, with n > 1. Note that we must have 
^2n+i < E2n-i and £^2n+2 < £^2n siucc larger n corre- 
sponds to successively better variational wave functions 
with successively lower energies. Now we assume for each 
rj that the corrections to the energies become smaller with 
increasing n such that 

^2n+l — ^2n+3 < ^2n " £^2n+2 < £^2n-l " £^2n+l (29) 

and we also assume that ^2n-i (^2n) is a monotonically 
decreasing (increasing) function of 77 around the transi- 
tion like in Fig. 6, i.e. near a given point r/o we can write 



0.75 



N 0.5 



0.25 




E2n-l{v) 



E2n{Vo) -\- pM{V-m) 
E2n-l{Vo) - PP{V-V0) 



FIG. 7: Residue of the attractive polaron for equal masses 
calculated from the variational wavefunctions ^5(0) (solid 
line) and ^3(0) (dashed). The dashed vertical lines repre- 
sent the polaron-molecule transitions calculated within the 
two approximations. 



where pp^pu > 0. Then for successively better approxi- 
mations for the polaron-molecule transition: 

^2n(r/+)=^2n-l(r/+) 
E2n{V~) = E2n+l{v~) 
^2n+2(r/0=^2n+l(r/0 

one can prove that r]~ < r]^ < r]^ . Therefore the P^-Mq 
transition will lie between the P^-M^ and P3-M4 transi- 
tions, and so on, thus proving that the results (28) pro- 
vide upper and lower bounds. 

Finally, we also calculate the wavefunction overlap 
with the free impurity, the residue Z. As illustrated in 
Fig. 7, the residue also agrees quite well between the 
two ansatze, even across the regime of strong many-body 
effects. Note that the residue is finite at the polaron- 
molecule transition, indicating that the residue jumps to 
zero discontinously when the impurity binds an extra t 
particle. This is a signature of a sharp "first-order" bind- 
ing transition like in 3D [19, 50]. To our knowledge, our 
result for IP5 ) is the best estimate for the polaran residue 
thus far presented in the literature. 



VII. FINITE DENSITY 

The high polarization limit of the phase diagram for 
the spin-imbalanced Fermi gas involves a finite density 
of impurities, and thus the question arises whether the 
single-impurity transitions are thermodynamically stable, 
i.e., whether they are preempted by first-order transitions 
in the thermodynamic limit. Typically, one requires the 
interactions between dressed impurities in order to as- 
sess this scenario; for instance, an attractive interaction 
between bosonic dressed impurities (e.g. dressed dimers 



or tetramers for a spin-down fermionic impurity atom) 
implies that the dressed impurities are unstable towards 
collapse into a region of higher density. However, the 
approach we will employ here is to determine whether 
the conditions for phase separation between the super- 
fluid (SF) and fully-polarized normal (N) phases are ever 
satisfied. To this end, we assume that the superfiuid is 
unpolarized and exploit the equal-mass equation of state 
for the BCS-BEC crossover derived from QMC calcula- 
tions [32]. This assumption is reasonable since no polar- 
ized superfiuid is observed in the mean-field calculations 
for the 2D Fermi gas [51]. 

The onset of phase separation in a Fermi gas near full 
polarization corresponds to the following conditions for 
the pressures and chemical potentials in each phase: 



PSF = Pn • 



= E. 



(30) 



where E is the energy of a single ^ impurity immersed in 
a spin-polarized Fermi gas. At the mean- field level, this 
is very easily carried out by setting E = 0^ i.e. assum- 
ing that the impurity atom is non- interacting, and then 
analyzing the minima of the thermodynamic potential 
r^(A,/icr) = —pV^ where A is the mean-field superfiuid 
order parameter. Identifying the points at which we have 
degenerate minima at A = and A 7^ corresponds 
exactly to satisfying the conditions (30). Using the 
mean-field approach, we obtain a first-order transition at 
Sb/Ef^ = Vl + r — 1. Below, we obtain a more accurate 
result for the case of equal masses = = m. 

Clearly, the pressure and spin-up chemical potential in 
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the fully-polarized normal phase are known exactly: 



Pn 



167rm 



E 



Ft 



2m 



(31) 
(32) 



Thus, we just need to estimate the binding energy E = 
AEp^ to determine fi^ . We obtain estimates from both 
wave functions \Ps) and IP5). 

The pressure and average chemical potential for the 
unpolarized superfluid can be written: 



PSF 



87rm 

SF 

SF 



Cp(^) 



2m 



(33) 
(34) 



where Cp{^)^ Cij,{^) interpolating functions that can 
be determined from the QMC equation of state [32] , and 
jy = ln(/ci?a2D) = 77 + ln(/ci?//ci?>^). The Fermi momentum 
kp of the superfluid is generally different from that in the 
fully-polarized phase. 

By equating pressures and chemical potentials in each 
phase, we obtain the coupled equations 



kp^ 
kp 



kp^ 
kp 



(35) 



which must be solved for 77 and kp^/kp- For the SF-N 
transition, we thus get within the two approximations for 
the impurity wave function 



IPs] 



rj 



-0.69, Sb/Ef^ 
-0.90, eB/Ept 



7.9 
12.2 



(36) 



By comparing these values with (28), we see that this 
first-order transition (and the concomitant phase separa- 
tion) preempts the single- impurity transition in the ther- 
modynamic limit. This implies that the single-impurity 
transition effectively corresponds to a spinodal line for 
the SF-N transition, and this could be the case for much 
of the polaron- molecule transition line in Fig. 5. Note, 
however, that the results (36) strictly apply to zero tem- 
perature and require the existence of a superfluid. Such 
a SF-N transition may be destroyed by strong thermal 
fluctuations and so we can envisage a scenario where we 
only have the single-impurity transition existing at flnite 
temperature. Future work is required to distinguish be- 
tween these possible scenarios. 



VIII. VARIATIONAL DESCRIPTION OF 
METASTABLE STATES 

The variational approach can also be generalized to 
study metastable excited states such as the repulsive po- 
laron [42, 43, 52, 53]. Rather than minimizing the energy 



of the variational wave function, one instead needs to con- 
struct its equations of motion by minimizing the "error" 
quantity [54, 55] 



for all allowed variations of the unknown function idtil^{t)^ 
where 5 = idt — H. For the polaron wave function (11), 
this yields the same equations (12) and (13), but with 
energy E replaced with idt. Now for a metastable polaron 
state with a long lifetime, we consider time-dependent 



amplitudes of the form 4^^(t) = a^^\0)e-'^^-^^ and 
a[p)(t) = Q;^^^(0)e-^^^+^^ where the decay rate F < 



IE"!. Note that we require to grow exponentially 
with time while Qq^^ decreases exponentially so that the 
normalization condition \a^\t)\'^ + Zlkq l^kq^ (^) P = ^ 
is preserved up to leading order in T/E. Inserting these 
amplitudes into the dynamical equations for the polaron 
then yields: 



E-iT-e 



Pi 



q 

^{p,E- 



k 

■iV) 



E, 



pqk 



iV 



which corresponds exactly to the condition for the quasi- 
particle pole within the diagrammatic approach [53]. 
Thus, we obtain the usual equations for the energy and 
the decay rate of the repulsive polaron: 

^ = ep; + 5R[S(p,E + iO)] 

r--^[s(p,E + io)] 

In principle, this approach can also be used to study 
the metastable states proximate to the binding transi- 
tions [50] . Note that energy and momentum conservation 
restricts the possible decay channels and means that a 
variational description which includes the metastability 
needs necessarily to include dressing by extra particle- 
hole pairs. 



IX. 



CONCLUDING REMARKS 



This work presents a thorough investigation of the 
highly polarized limit of a 2D Fermi gas, which can be 
modelled as a J. impurity in a t Fermi sea. We have an- 
alyzed the possible states (polaron, molecule, trimer...) 
that the impurity can form and we have constructed the 
single-impurity phase diagram using two different levels 
of approximation for the impurity wave function. The 
simple "undressed" wave functions in Sec. Ill give us in- 
sight into the few-body correlations of the system and 
allow us to make contact with standard mean-field ap- 
proaches for the spin-imbalanced Fermi gas. The wave 
functions dressed with one particle-hole excitation in 
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Sec. IV contain the correct first-order correction to the 
impurity energy due to interactions with the medium, 
and they yield a phase diagram that should be qualita- 
tively, if not quantitatively, accurate across the full range 
of interactions 77. Indeed, our study of the polaron state 
with two particle- hole pairs in Sec. VI suggests that our 
variational approach is reasonable even in the regime of 
strong interactions 77 ~ 0. To better model experiment, 
the variational calculation may also be extended to prop- 
erly include the transverse harmonic confinement — for 
further details, we refer the reader to Ref. [24]. 

Our results show that the trimer is strongly favored 
by the Fermi sea since the FFLO molecule is unstable to 
binding an extra t fermion, like in 3D. However, there 
is still the possibility of observing the FFLO phase at 
lower mass ratios 1.8<r<2.5. A remaining ques- 
tion is what happens to the trimer phase when there is 
a finite density of impurities. Naively, one might expect 
a Fermi liquid of trimers if they are sufficiently tightly 
bound, but the fact that they have finite angular mo- 
mentum could impact the properties of this liquid phase. 
There is also the question of whether the single- impurity 
binding transitions are thermodynamically stable — we 
have already shown that the polaron-molecule transition 
for equal masses is preempted by a SF-N transition at 
zero temperature. Such single-impurity transitions may 
nonetheless survive at finite temperature, in which case a 
finite density of molecules would correspond to a normal 
phase of preformed pairs. 

Finally, we have computed rigorous upper and lower 
bounds for the position of the polaron-molecule transi- 
tion in the case of equal masses. This should provide a 
useful benchmark for future experimental and theoretical 
work on this topic. 
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Appendix A: Bound states in the N -\- 1 problem 

In this Appendix we demonstrate how the equations for 
the binding energy of a composite consisting of N ma- 
jority particles and 1 minority particle may be obtained 
by a diagrammatic technique. This is an alternative to 
the variational method described in the main text. 

Consider first the undressed wavef unctions, the subject 
of Section III. For simplicity of notation, we will here 
consider the quasiparticle state at rest. The generaliza- 
tion to finite momentum may be easily achieved. The 
t particles have momenta ki, . . . , and corresponding 
energies ^kit? • • • ^^kivt' while the impurity has momen- 
tum — hi and energy E — ^k^t' such that the total 
energy is simply E. 

The sum of diagrams with N -\- 1 incoming particles 
where the impurity interacts first with the t particle of 
momentum ki will be denoted /ks-.-k^v- The function / 
does not depend on ki since the initial interaction de- 
pends solely on the total momentum of the two parti- 
cles. As the majority particles are fermions, / is anti- 
symmetric in its indices. The specific ordering of indices 
is of course arbitrary, however once an ordering is chosen 
it must be kept throughout the calculation as this cor- 
responds to a choice of ordering of operators in Wick's 
theorem. 

The occurence of an A/" -h 1 particle bound state with 
binding energy E corresponds to a singularity of / at this 
energy. For the polaron, / is simply the bare impurity 
propagator and the energy is that of the non-interacting 
impurity. For the molecule, / is the pair propagator in 
the medium, denoted T2. For N > 2 the singularity 
appears from the summation of an infinite number of di- 
agrams, and may be found by solving the integral equa- 
tion illustrated in Fig. 8: The initial interaction between 
the impurity and the particle with momentum ki is de- 
scribed by a pair propagator. Subsequently, the impurity 
interacts with another of the initial particles. Thus the 
right hand side contains N — 1 terms and the function / 
satisfies the integral equation 



/k2 



1 

^ki...kiv ~ ^0 



= E 

ki 



/kik3...kjv ~l~ /k2kik4...kAr + • • • + /k2...kAr_iki 



^ki. 



iO 



(Al) 



where ^ki...kjv = -E -\- eki+...+kA,^ + ^k^t 
factor — iO acts to slightly shift the energy pole into the 
lower half of the complex plane. The quantity in brackets 
on the left hand side is the inverse pair propagator. 



Appendix B: The dressed N -\- 1 problem 

We turn now to the states dressed by one particle- 
hole excitation as investigated in Section IV. Again, we 
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-ki - . . . - kiv 



k'l 



FIG. 8: The diagrams which lead to the binding energy of the N -\- 1 composite. Black dots on vertices indicate that the 
corresponding particles interact first inside the vertex /. 
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FIG. 9: The equations which give the binding energy of the N composite including dressing by one particle-hole pair. Black 
dots on vertices indicate that the corresponding particles interact first inside the vertex, (a) The equation for the vertex J, 
there are N — 1 diagrams of the first type on the r.h.s. corresponding to interactions between the impurity and one of the 
particles with momentum k2, . . . , Icat and 1 of the second type, (b) The equation for G. The first two types of diagrams on the 
r.h.s. appear N times corresponding to interactions between the impurity and one of the particles with momentum ki, . . . , kiv. 
Dashed lines indicate that these are the loops which are closed upon insertion of G in the diagrams in (a). 



construct first the sum of all diagrams with incoming 
t particles in which the impurity interacts first with the 
t particle of momentum ki. The kinematics is chosen 
as above. Now we allow the vertex to be dressed by one 
particle hole pair and this new vertex is denoted J\c2...\cn- 
Again, the initial interaction is described through a pair 
propagator. Then, in addition to terms of the same form 
as above where the impurity interacts next with another 



of the initial particles, it may also interact next with a 
particle from the Fermi sea. This is illustrated in Fig. 9a 
where the sign arises from the fermion loop. The vertex 
Gqki...kiv is of diagrams with one incoming 

I particle and N -\- 1 incoming t particles, the initial 
interaction being between the impurity and the particle 
having momentum q and energy ^^-f- The equation in 
Fig. 9a is 



ki ^1 



iO 



E 

ki 



ik3...kAr ik4...kAr + . . . + Jk,... KiV-lKi 



G. 



qki...krt 



ki...kA 



iO 



kiq 



iO 



(Bl) 



The vertex G is constructed as follows: The initial inter- mentum q is again described by a pair propagator. Af- 
action between the impurity and the particle with mo- 
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ter this interaction, the t particle which partakes in this 
repeated interaction is either disconnected or connected 
from the remaining particles. If it is connected then the 
impurity interacts next with one of the particles with mo- 
mentum ki, . . . , k^v, as the vertex is only dressed by one 
particle hole pair. This interaction is again described by 



the vertex G. If disconnected, the impurity can interact 
first with one of the particles with momentum ki , . . . , kAr 
and this interaction is given by the vertex J, or it can 
interact first with a particle from the Fermi sea. This 
results in the equation for G: 



qki 



qki...kjvk 



iO 



_ '^-^ Gqkk2...kAr ~^ Gqkikks.-.kjv H" • • • H" Gqki ...kAr_ik 
~ V ^qki...k^k-iO 



«/k2...kAr — «/kik3...kjv + -^kik2k4...kAr — • • • + ( — 1)^ ^ «/ki ...kjv-i Sq' Gq^ki. 



^ki...kiv ~ ^0 



(B2) 



where ^qki...kivk = -E + eq_ki...-kiv-k; + Z^i ^^k,t + 
^kt ~ ^qt- (-^^) illustrated in Fig. 9b where the 
dotted lines, indicating how the loops are closed, should 



only be taken as a guide to the eye, as the integration 
over momentum in these loops is not performed until the 
insertion of this equation in Eq. (Bl). 
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